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Abstract 



A search for the decay K$ — > 7r°e + e _ has been made by the NA48/1 experiment 
at the CERN SPS accelerator. Using data collected during 89 days in 2002 with a 
high-intensity K$ beam, 7 events were found with a background of 0.15 events. The 
branching fraction BR(K S -^ir°e + e-, m ee > 0.165 GeV/c 2 ) = (3.0+^(stat) ± 
0.2(syst)) x 10 -9 has been measured. Using a vector matrix element and a form 
factor equal to one, the measurement gives BR{K$ — > 7r°e + e~) = (5.8^24) x 10 -9 . 



1 Introduction 

When not forbidden by CP-conservation, the K — > iree decay can proceed via single 
photon exchange. This is the case for K$ and charged kaons, while the Kl decay - barring 
a small CP-conserving contribution - is CP- violating. 

The rate of K s — > 7r°e + e~induced by the electromagnetic interaction was predicted 
in ref. [T] to be BR(K S -> 7r°e + e-) = 5.5 x 1(T 9 . 

The theoretical aspects of the decay Ks — > ir°e + e~ were studied to leading order in 
the chiral expansion in ref. [2J dj and the implications of this decay with respect to the 
search for CP-violation in rare kaon decays were investigated in ref. [3] and re-examined in 
ref. Further study beyond leading order was presented in ref. [Hj, where the branching 
fraction for K s — > n°e + e~ was expressed as a function of one parameter a s : 

BR(K S -> 7r°e + e-) = 5.2 x 1(T 9 a 2 s . (1) 

For Kl — > TT°e + e~ , CP- violating contributions can originate from: 

a) K° — K° mixing via a decay of the CP-even component of the Kl (Ki) into n°e + e~ . 
This indirect CP-violating contribution is related to the Ks branching ratio: 

BR(K L AV) = T ^BR(K S A+e-) K gg^£_^£) f2) 

Ts OOU 

b) direct CP-violating contribution from short distance physics via loops sensitive to 
Im(X t ) = Im(V td V t * s ). 

The indirect and direct CP-violating contributions can interfere and the expression 
for the total CP- violating branching ratio of Kl —>■ 7r°e + e~can be written as [6]: 

BR(K L -> 7r°e + e-)cPv x 10 12 ~ 15.3 a| - 6.8 a s (^^) + 2.8 (^^) ' (3) 

As shown in eq. EJ the sensitivity to Im(X t ) can also come from the interference 
term depending on the value of as- The theoretical predictions for Ks — > 7r°e + e~ do not 
provide firm constraints on Im(X t ) and a measurement or a stringent upper limit on as is 
necessary to progress further in the understanding of CP-violation in the Kl — > n°e + e~ 
decay. 

Currently, the upper limit of the BR(Kl — > 7r°e + e~) is 5.1 x 10~ 10 [7j. This to- 
gether with the present upper limit BR,(Ks — > 7i°e + e~) < 1.4 x 1CU 7 |B| gives a bound on 
Im(Xt) 05 but not competitive with respect to other constraints obtained from b-physics. 
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2 Data-taking 

2.1 Beam 

The experiment was performed at the CERN SPS accelerator, and used a 400 GeV/ c 
proton beam impinging on a Be target to produce a neutral beam. The spill length was 
4.8 s out of a 16.2 s cycle time. The proton intensity was fairly constant during the spill 
with a mean of 5 x 10 10 particles per pulse. 

Fig. ^ shows the modifications with respect to the previous K$ beam line described 
in ^Hj • The Kl beam line was blocked and an additional sweeping magnet was installed 
to cover the defining section of the Kg collimator. To reduce the number of photons in 
the neutral beam, primarily from 7r° decays, a platinum absorber 24 mm thick was placed 
in the beam between the target and a sweeping magnet, which deflected charged particles. 
A 5.1 m thick collimator, the axis of which formed an angle of 4.2 mrad to the proton 
beam direction, selected a beam of neutral long-lived particles (Kg, Kl, A , H°, n and 
7). On average 2 x 10 5 K$ per spill decayed in the fiducial volume downstream of the 
collimator with a mean energy of 120 GeV. 

2.2 Detector 

The detector was designed for the measurement of Re(e' / e) jTOJ • In order to minimize 
the interactions of the neutral beam with air, the collimator was immediately followed by 
a ~ 90 m long evacuated tank which was terminated by a 0.3% Xq thick Kevlar window. 
The detector was located downstream of this tank. 

2.2.1 Tracking 

The detector included a spectrometer housed in a helium gas volume with two drift 
chambers before and two after a dipole magnet with a horizontal transverse momentum 
kick of 265 MeV/c. Each chamber had four views (x, y, u, v), each of which had two 
sense wire planes. The resulting space points were typically reconstructed with a resolu- 
tion of ~ 150 fim in each projection. The spectrometer momentum resolution could be 
parameterized as: 

a p /p = 0.48% © 0.015% x p 

where p is in GeV/c. This gave a resolution of 3 MeV/c 2 when reconstructing the kaon 
mass in Kg — > 7i + n~ decays. The track time resolution was ~ 1.4 ns. 

2.2.2 Electromagnetic Calorimetry 

The detection and measurement of the electromagnetic showers were achieved with 
a liquid krypton calorimeter (LKr), 27 radiation lengths deep, with a ~ 2 cm x 2 cm cell 
cross-section. 

The energy resolution, expressing E in GeV, may be parameterized as [TTj : 

a{E)/E = 3.2%/VE © 9%/E © 0.42% 

The transverse position resolution for a single photon of energy larger than 20 GeV 
was better than 1.3 mm, and the corresponding mass resolution at the ir° mass was ~ 
1 MeV/c 2 . The time resolution of the calorimeter for a single shower was better than 
~ 300 ps. 
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2.2.3 Scintillator Detectors 

A scintillator hodoscope was located between the spectrometer and the calorimeter. 
It consisted of two planes, segmented in horizontal and vertical strips and arranged in four 
quadrants. Further downstream there was an iron-scintillator sandwich hadron calorime- 
ter, followed by muon counters consisting of three planes of scintillator, each shielded by 
an iron wall. The fiducial volume of the experiment was principally determined by the 
LKr calorimeter acceptance, together with seven rings of scintillation counters used to 
veto activity outside this region. 

2.2.4 Trigger and Readout 

The detector was sampled every 25 ns with no dead time and the samples were 
recorded in a time window of 200 ns encompassing the event trigger time. This allowed 
the rate of accidental activity to be investigated in appropriate time sidebands. 

The event trigger for the signal K$ — » 7r°e + e~ had both hardware and software 

parts: 

— The hardware trigger [T2] selected events satisfying the following conditions: 

hit multiplicity in the first drift chamber compatible with one or more tracks; 

hadron calorimeter energy less than 15 GeV; 

electromagnetic calorimeter energy greater than 30 GeV; 

the centre of energy of the electromagnetic clusters (see eq. El below) less than 

15 cm from the beam axis; 

the decay occurring within six K$ lifetimes from the end of the collimator; 
no hits in the two ring scintillator counters farthest downstream. 

— The software trigger required: 

at least two tracks in the drift chambers and two extra, well-separated clusters 
each with energy greater than 2 GeV; 

the tracks projected from the drift chamber, after the magnet, had to match 
to clusters in the LKr within 5 cm; 

the tracks had to be compatible with being electrons or positrons using the 
condition that the ratio E/p, between the cluster energy in the LKr, E, and 
the momentum measured with the drift chambers, p, had to be greater than 
0.85. 

a cluster separation of more than 5 cm was required to limit the degradation 
of the energy resolution due to energy sharing between closely spaced clusters. 
The events that satisfied the trigger conditions were recorded and reprocessed with 
improved calibrations to obtain the final data sample. 

2.3 Event selection 

For the analysis of the data, signal and control regions were defined. These regions 
were masked while the cuts to reject the background were tuned using both data and 
Monte Carlo simulation. 

The signal channel Ks — > 7i°e + e~ required the identification of an electron and a 
positron accompanied by two additional clusters in the LKr. 

Tracks reconstructed from the spectrometer which matched an LKr cluster were 
labelled as an electron or positron by requiring three conditions to be met: no more than 
3 ns difference between track time and cluster time; 0.95 < E/p < 1.05; and less than 
2 cm between the projected track and the cluster coordinates in the LKr. 
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We define At to be the difference between the average time of the two clusters 
associated with tracks and the average time of the two neutral clusters. Events were 
accepted if \At\ < 3 ns. 

Events with extra tracks or extra clusters within 3 ns of the average time of the 
tracks or clusters and with an energy larger than 1.5 GeV were rejected. To minimize the 
effect of energy sharing on cluster reconstruction, a minimum cluster separation of 10 cm 
was imposed. In addition, a distance greater than 2 cm between the impact points of the 
two tracks at the first drift chamber was required. 

Four quantities related to the decay vertex were computed. 

— neutral vertex. 

The vertex position was computed from the energies and positions of the four clus- 
ters in the LKr according to 

Zneutral = %LKr ~ \J ^i,j>iEiEjdfj/ M K (4) 

where zlkt is the longitudinal position of the front face of the LKr; M K is the kaon 
mass, Eij is the energy of the {i,j) th cluster and dij is the distance between clusters 
% and j. In the case of the photons these are the x, y shower positions in the LKr. 
For the e ± tracks, in order to cancel the deflection due to the dipole magnet, the 
(x,y) positions were calculated by extrapolating the tracks from their positions in 
the first two drift chambers to the face of the LKr. The x and y coordinates of the 
neutral vertex were found by extrapolating the position of each track before the 
magnet to the position of z ncutral . The average of the two measurements was taken 
as the (x, y) vertex position. 

The neutral vertex was used to compute the invariant mass of the two photons, 

— charged vertex. 

The position of the charged vertex can be calculated using the constraint that the 
kaon decay should lie on the straight line joining the target and the point defined 

as (x CO g, y C og)- 



%cog 



EiXi)/Y, E i Vco 9 = E E lVl )/J2 Ei (5) 



where Ei, Xi and yi are the energy and positions of the i-th cluster. 
For each track, the closest distance of approach between this line and the track 
was found, giving two measurements which were then averaged to give the charged 
vertex position. 

The charged vertex was then used to compute m eeTY , the invariant mass of the four 
decay products. 

- 7T° vertex. 

The 7T° vertex position along the beam direction was computed in a similar way 
to the neutral vertex, but using only the two photon clusters and imposing the n° 
mass, M n o, instead of the kaon mass. 

— track vertex. 

The track vertex is at the position of the closest distance of approach of the two 
tracks. 

The z position of the tc° and track vertices had to be greater than 50 cm (one 
standard deviation) beyond the collimator exit in order to reject any interactions occurring 
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in the collimator. Assuming the observed event to be a kaon decay, the proper lifetime was 
computed from the position of the neutral vertex, taking the end of the final collimator as 
the origin. A cut at 2.5 Ks lifetimes was then applied. The kaon momentum was required 
to be between 40 and 240 GeV/c. 



3 Signal and Control regions 

The signal region was defined as : 

- |m 77 - M w o\ < 2.5 x a mil 

- |m ee77 - M K \ < 2.5 x a mee77 

To evaluate the resolutions, o" mee77 and <r m77 , we studied the channel K s — > ir ^ 1 *, 
for which we measured cr mee777 = 6.5 MeV/ c 2 and a m77 = 1 MeV/ c 2 respectively. These 
values were found to be in agreement with a Monte Carlo simulation based on GEANT 
[T3] . For the decay Ks — > 7r°e + e~, the Monte Carlo prediction of <J mee was 4.6 MeV/ c 2 
and this value was used in defining the signal region. The better resolution is due to the 
fact that the e + e~ opening angle is on average larger than for the decay Ks — > tt ^. 

The m 77 resolution, cr m77 , at the n° mass was found to be 1 MeV/ c 2 in agreement 
with the Monte Carlo simulation. 

A control region was also defined as: 

- 3 x er m77 < |m 77 - M T o \ < 6 x <j m77 

- 3 x a mee _ n < |m ee77 - M K \ < 6 x a m ^ n 

Both the signal and the control regions were kept masked while cuts to reject the 
background were studied. 



4 Background Rejection 

A large number of possible background channels was studied. These channels were 
of two types: 

— a single kaon or hyperon decay which reproduced an event falling into the Ks — > 7r°e + e 
signal region 

— fragments from two primary decays which happen to coincide in time and space 
and fall into the signal box. 

The background contribution from the channels considered was reduced by imposing 
additional requirements. 

A background source is from the decay Ks — > 7r 7r° where two photons from different 
7r°'s converted either internally (i.e. Ks — > ^d^d) or externally and one electron and one 
positron from different 7r°'s were outside the detector acceptance. In order to reject events 
from this source the invariant masses of the two electron-photon pairs, m e + 7l ,m e - 72 and 
m e + 72 ,m e - 71 , were computed using the charged vertex position. A priori, the combination 
corresponding to electron-photon from the same it has an invariant mass smaller than 
M n o. Thus events were rejected if both m n masses were measured to be smaller than 
M w o + 5. The constant 5 was chosen equal to 30 MeV/c 2 , which corresponded to ~ 10 cr me7 . 

Another source of background was due to 7r°7r° decays where one or more photons 
from a single tt° decay converted (either internally or externally). These decays are kine- 
matically constrained to have m e + e - < M n o and in order to reject this background the 
analysis was restricted to the event sample with invariant mass m e + e - > M n o + e. To de- 
termine e, we analysed the m ee distribution from data and compared it to a Monte Carlo 
simulation where the different components were identified. In fig. |21a we show the m ee 

l > tt° d is the Dalitz decay n° — > e + e~7. 
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distribution for data (full dots) and superimposed the contributions from of all relevant 
background sources. Above the ir° mass the tail of the m ee distribution falls rapidly to 
zero. The constant e was also chosen equal to 30 MeV/ c 2 and the analysis was therefore 
restricted to the region m e + e - > M n o + 30 MeV/c 2 = 165 MeV/c 2 , where 7 conversions 
or decays from a single 7r° give a negligible contribution to the background. This was con- 
firmed from the analysis of events with same sign-tracks. This sample contained events 
where both photons from a single ir° converted and both the electrons or the positrons 
were in the acceptance. The m e ± e ± distribution is shown in fig. Elb, where data and Monte 
Carlo are compared. No events with m e ± e ± > 165 MeV/c 2 were found. 

To reject the background due to electron bremsstrahlung, the invariant mass of any 
cy combination was required to be larger than 20 MeV/ c 2 . 

The background from S° — > Air and A — > pir~ decays was reduced to a negligible 
level by exploiting the large momentum asymmetry in both the Att° and the pn~ final 
states. Ks — > 7r°e + e~ candidates were required to have (P^ — P n o)/ (P^ + P n o) smaller than 
0.4 or (P p — P T -)/(P p + P n ~) smaller than 0.5. A similar cut was used to remove H° and 
A. 

The possibility of proton and pion misidentification as e ± was considered and final 
states which contained these particles were found to make a negligible contribution to the 
background after the application of the E/p requirement. 

5 Estimate of the Residual Background 

After the selection outlined above three sources of background were found to be 
non-negligible: 

1. K L ,s -> e+e"77 

The e + e~77 component was measured using Kl data from the 2001 run, in which 
the number of Kl — > e + e~77 decays was ~ 10 times the sum of the Kl and 
K s — > e + e~77 expected in the present experiment. The distribution of m ee77 versus 
m 77 for these events is shown in fig. El Using a linear extrapolation from the low 
m 77 region to the signal region, the background from this channel was estimated to 
be 0.08±^ events. 

2. K s - 7r°,7r° 

This was evaluated using full Monte Carlo simulation for a sample which was 30 
times greater than the data, and the background was estimated to be less than 0.01 
events in the signal region. 

3. Accidental backgrounds. 

This component was studied using data with the timing requirements relaxed. 
Events in the time sidebands, satisfying all the other cuts, were used to extrap- 
olate the background from the control to the signal region. A further correction was 
applied to account for the background shape in the m 77 versus m e + e - 77 plane as 
predicted by a simulation. 

The contribution due to this component was 0.071qq3 events in the signal region. 

Other sources of background were considered, for instance that due to resonances 
produced by a single proton in the target, and decaying to a pair of kaons or a KA pair 
in the fiducial region. These contributions were found to be negligible. 

With all the cuts applied, the control region was unmasked to estimate the final 
background contribution to the signal. No events were found in the control region, con- 
sistent with the background prediction of 0.33 events. Only one background event was 
found in a much larger region (corresponding to 17 x a mK and 20 x a m ). The background 
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estimate is summarised in Tabled 

The resulting estimate of the total background in the signal region was 0.15^q;q4 
events. 

6 Normalisation. 

The trigger efficiency was measured using a control sample of ~ 8.6 x 10 6 K$ — > 7r°7r^, 
decays, which differed topologically from K$ — > 7r°e + e~ only in having an extra photon. 
This sample was collected with the same trigger chain. The trigger efficiency, measured 
with a sample of triggers collected requiring minimal bias conditions, was found to be 
99.0%. The acceptance, including the selection criteria, was found to be 3.3% for 7r°7r£,, 
evaluated using a Monte Carlo simulation. The Monte Carlo simulation was found to be in 
good agreement with the Ks — > 7r°7r£, data. To obtain the Ks —>■ 7c°e + e~ branching ratio, 
the Ks flux was calculated using the channel K$ — > tt ^ for normalisation, which was se- 
lected using the same trigger. Using the value for the branching ratio BR(Ks — > tt ^) = 
7.43 x 10 -3 Jl] the Ks flux was calculated, for kaon energies between 40 and 240 GeV/ c 
and kaon lifetimes between zero and 2.5 Kg mean lifetimes from the collimator exit. The 
total number of Ks decaying within the fiducial volume was (3.51 ± 0.17) x 10 10 . 

7 Result 

When the signal region was unmasked seven events were found (fig. ^). With an 
expected background of 0.15^o;q4 events, this corresponds to a sig nal of 6.85±?;|. The 
probability that such a signal is consistent with background is ~ 10~ 10 . We therefore 
interpret the signal as the first observation of the Ks — > 7i°e + e~ decays. 

Fig. [5] shows the m 77 and the m e + e - 77 distributions of the events compared to the 
detector mass resolutions. In Table El some of the kinematical quantities for each event 
are summarised. 

In order to calculate the Ks — > n°e + e~ acceptance, the amplitude for the decay was 
needed. This was taken from the Chiral Perturbation Theory prediction given in [E] , and 
is of the form: 



A[K(k) - n(p)e + (p + )e-(p_)] = W{z)(k + P yu l ( P -) lfl v l {P+) (6) 

rriK yttr) 

where k, p, p + and p_ are the four-momenta of the kaon, pion, positron and electron 
respectively; is the kaon mass; W(z) is the electro- magnetic transition form factor, 
with z = (k — p) 2 /mK 2 - As a consequence of gauge invariance, the form factor dependence 
on z vanishes to lowest order and therefore can be represented as a polynomial. For Ks 
decays, the form factor W(z) was approximated to W(z) ~ a + b x z [HJ. 

The a and b parameters have recently been measured for charged kaons, and the 
ratio a/b found to be 1.12 [T5j . 

The m ee distributions resulting from W(z) = 1 and W(z) = a + b x z are shown in 
fig. 0a. 

The overall Ks — > 7r°e + e~acceptance depends on the form factor. To remove this 
form factor dependence, an acceptance was calculated for each event using fig. |3b, where 
the acceptance is given as a function of m ee , using W(z) = 1. The values are given in 
Table |2j The average geometrical acceptance of the 7 events is 0.15, while the average 
analysis efficiency is 0.44, which results in an average efficiency of 0.066 ± 0.004. 
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From the Kg — > tv tt% flux and the signal of 6.85 events, the branching ratio for 
m ee > 0.165 GeV/c 2 was computed: 

BR(K S 7r°e + e-, m ee > 0.165 GeV/c 2 ) = (3.0l};ij(stat) ± 0.2(syst)) x 10~ 9 . 

The quoted uncertainties correspond to a 68.27% confidence level ^E]- The system- 
atic uncertainty includes the uncertainty of the flux measurement and of the acceptance. 

8 Discussion 

In Chiral Perturbation theory the BR(Ks — > 7r°e + e~) is related to the parameter 
as, which measures the strength of the indirect CP-violating term in K~l — > 7r°e + e~decay 
as explained in [Hj and eq. 

Using a vector matrix element with no form factor dependence, the measured 
branching ratio was extrapolated to the full m ee spectrum to obtain: 

BR(K S -> vr e + e-) = (5.8^|(stat) ± 0.8(syst)) x 10~ 9 

The systematic error is dominated by the uncertainty in the extrapolation due to 
the form factor dependence. 

It was then possible to extract the parameter \as\- 

\a s \ = (1.06lg|f(stat) ± 0.07(syst)) 

The measurement of as allows the branching ratio BR,(K~l — > 7r°e + e~) to be pre- 
dicted as a function of Im(X t ) to within a sign ambiguity (see eq. EJ). The effect of the 
sign ambiguity can be seen in fig. 0a. 

Alternatively, as shown in fig. 0b, by using the global fit value for Im(X t ) = (1.30 ± 
0.12) x 10~ 4 obtained from b-decay PH, BR( J K" L — > 7r°e + e ) can be expressed as function 
of \as\- 

Using the measured value of \as\ and the global fit for Im(X t ), eq. El reduces to: 

BR{K L -> 7r°e + e-) C pv ^ (17.2 indircct ± 9.4 intcrfcrence + 4.7 direct ) x 10~ 12 . 

The CP conserving component can be obtained from the study of the — > 7r°77 
decay. A measurement made by the KTeV collaboration [THj found BR(K L — > 7i°e + e^)cpc — 
(1 — 2) x 10 -12 . A more recent measurement quoted BR(K L — ► 7c°e + e~)cpc = 0.47^g; 2 g x 
10~ 12 JH] suggesting that the CP-conserving component is negligible. 

Given the measured value of as the direct CP violated component predicted from 
the Standard Model is small with respect to the indirect component. If the sign of as turns 
out to be negative then BR(i^^ ~~ * ir°e + e~) retains some sensitivity to lm(X t ) through the 
interference term. 
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Source 


control region 


signal region 


K s -> n° D ir° D 


0.03 


< 0.01 


K L ee77 


0.11 


0.08 


Accidentals 


0.19 


0.07 


Total background 


0.33 


0.15 



Table 1: Summary of the background estimate. 



Event no. 


Ks momentum 


t/t s 


m ee 


Acceptance 




GeV/c 




GcV/c 2 




1 


84.6 


0.74 


0.291 


0.058 


2 


128.2 


0.50 


0.267 


0.066 


3 


114.1 


1.02 


0.173 


0.084 


4 


83.9 


2.09 


0.272 


0.066 


5 


130.8 


1.46 


0.303 


0.052 


6 


121.2 


1.49 


0.298 


0.058 


7 


94.2 


1.64 


0.253 


0.075 



Table 2: Kinematical quantities of the seven events found in the signal region. 




Figure 1: View of the 2002 modifications to the beam line. 
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• 2002 data 

I 7t°4(m ee >0.09GeV) 

D 7t°7t° + conversions 
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D 7t°Jt|?D 



: U 7t Jt DD | | |j 
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m ee (GeV/c 2 ) 




Figure 2: Distributions of m ee after all the cuts have been applied. Superimposed we show 
the Monte Carlo predictions from all important sources. Figure a) shows the components 
with opposite-sign tracks; figure b) shows the ones with same-sign tracks. 



0.65 




(GeV/c") 



Figure 3: Scatter plot of m eeil versus m 77 for events selected as Kl — > e + e 77 in the 
2001 data. The boxes are representing the 3a and 6a regions. 
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Figure 4: Scatter plot of m eery versus m ee (a) and m eery versus m 77 (b) for events passing 
all the cuts described in the text. The regions of 3a and 6a are shown. 
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Figure 5: m 77 (a) and m ee77 distributions (b) for the 7 events found in the signal region. 
The expected Gaussian mass resolutions are superimposed (solid line). 
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m ee (GeV/c 2 ) m ee (GeV/c 2 ) 



Figure 6: m ee distributions from Monte Carlo, with and without the form factor (a); 
acceptance as a function of m ee (b). 




Figure 7: Branching fraction of Kl — > 7t°e + e as a function of Im(\ t ) (a), and as a function 
of the parameter as (b). 
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